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I. INTRODUCTION 

The static properties of mean-field spin-glasses have been investigated in the past mainly by means of two different 
approaches: standard thermodynamics, which through the replica method jl], ||, || [|, ||] aims to compute the partition 
function and the equilibrium free energy density of the system, and the Thouless-Anderson-Palmer (TAP) method 
H @i |[ H 0' Ull' which introduces a mean- field free energy Ftap, function of the set m of local magnetizations, 
m = {rrii^^Z^. The local minima of FTAp(m) are identified with the metastable states of the system, which at the 
mean-field level are well defined. At low temperatures the total number Af of minima of the TAP free energy becomes 
exponentially large [ffl p[ [Tof, and the density of states is given by, 

P(f) ~ . (1) 

In this expression £(/) is the complexity of the TAP states with free energy density between / and / + df. 

The existence of two apparently different methods for the study of the static properties of mean-field spin-glasses 
clearly poses a problem of consistency. Both these methods are in principle correct, and the equivalence between the 
two approaches should hold in general. A direct argument in its favour is the cavity approach (pj, |l4| : on the one hand 
the cavity approach is known to be mathematically equivalent to the replica approach; on the other hand the first 
step of the cavity approach consists in clustering the configurations into states which are nothing but the TAP states. 
From the cavity analysis, one can expect a priori the existence of a formal relationship between the replica results 
computed at a given order k of replica symmetry breaking, and the results of TAP solutions, correctly weighted and 
with an average computed at the order k — 1. Unfortunately, when considering the detailed results published so far, 
there seemed to exist a discrepancy between the two approaches. For instance the number of TAP solutions computed 
directly at the order k = (annealed average) in |Q, disagrees with the one obtained indirectly through the replica 
method with one step of replica symmetry breaking (i.e. a k = 1 computation) p"4j . This apparent discrepancy has 
led us to give a careful look at this consistency problem. We show here that the TAP approach must be reconsidered: 
using a supersymmetry, one finds that the two approaches are in fact fully consistent. 

From a general point of view, the problem of consistency between the TAP approach and standard thermody- 
namics can be posed in a twofold way. First, the two methods must agree in the calculation of the thermodynamic 
quantities that can be computed in both frameworks. This requirement is a strong form of consistency. At a weaker 
level, regarding objects which are inherently defined only within the TAP method, we may still expect to find some 
consistency, and even some formal connections, with the standard thermodynamic approach. 

The strong consistency basically asks that the equilibrium free energy computed within the two approaches must 
be the same. Given that Ftap is a function of the local magnetization, it is not obvious how to compute the partition 
function in this framework. More specifically, due to the large degeneracy of TAP minima, there is the problem of 
how to weight correctly these states. This problem was solved by De Dominicis and Young |l2[ (DDY), who proposed 



2 



to weight each TAP state a with, 

w a = e X p[-{3F TAP (m a )} / Z TAP . (2) 

The strong equivalence of the TAP approach to the standard thermodynamic approach is therefore encoded in the 
equality (which should hold in the thermodynamic limit) logZ = -k log Ztap, that is, 

i log exp[-/3ff(a)]^j = 1 log ^ exph/3F TA p(m")]^ . (3) 

Proving relation (^) in general is highly nontrivial, especially because in spin-glasses we can only compare averages 
over the disorder, and thus, from a practical point of view, we need to compare the average of the replicated partition 
functions, which will typically require some tricky integration over complicated order parameters. This program has 
been carried out by DDY for the Sherrington-Kirkpatrick (SK) model Q). Their conclusion is that relation (||) is 
indeed verified, provided that the complexity of the dominant TAP states at any temperature is zero. This restrictive 
hypothesis is satisfied in the SK model, but not in other models where the TAP approach is nevertheless used. Thus, 
it may seem that either condition (||) is not broadly valid, or that a general proof of it must give up the hypothesis 
used by DDY. 

In the present paper we show that the proof of (0) given by DDY for the SK model is valid with no extra hypothesis. 
The key point is that the TAP partition function can be written in an integral form, with an action which is invariant 
under the Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry |l5|, This symmetry provides the mathematical 
relations needed to prove equation (||), without the need to invoke any further hypothesis on the complexity of the 
states. 

Concerning the weak consistency of the two methods, the key object is the complexity £(/). Even though there is no 
obvious way to compute £ in the standard thermodynamic framework, we may still impose an important consistency 
condition: at low temperatures (below the static transition) we expect the equilibrium thermodynamic states to have 
the same free energy density as the lowest TAP states. In other words, we require that at low enough temperatures 
there must be no static contribution of the metastable states. This condition is encoded in the equation, 

fe g = h , (4) 

where f eq is the equilibrium free energy density computed in the standard thermodynamic framework, and /o is 
defined by the TAP relation £(/o) = 0. Apart from this minimal requirement, one may ask whether in the TAP 
calculation of £ emerge some deeper formal connections with the standard calculation of the free energy. This fact, of 
course, seems more than likely if the strong consistency condition (0) holds, even though it is in general not obvious. 
In the p-spin spherical model ]l7j , relation (jl]) is satisfied jl8|, [l9|, Moreover, for this same model it has been 

shown in p2] that by means of the BRST supersymmetry the saddle point equations involved in the TAP calculation 
of £ at / = fo become identical to the static equations in the replica approach. In the SK model, however, relation 
(|J) is not easy to check Q [§, 0, [HJ O, ^3) , since it requires a full replica symmetry breaking (RSB) calculation of the 
complexity. Moreover, up to now there was no formal connection between the TAP complexity and equilibrium free 
energy in the SK model. 

The fact that the BRST supersymmetry plays a crucial role in the p-spin, suggests that even in the SK model this 
symmetry may help to prove (|j), and discover other formal connections between the two methods. This is what we 
show in Section IV, where we revisit the annealed calculation of the SK complexity of Bray and Moore [Q (BM), 
and thanks to the BRST supersymmetry show that at / = fo this calculation becomes equal to the standard static 
calculation of the free energy, at one step level of replica symmetry breaking. This result comes hardly as a surprise, 
once the strong consistency (||) of the two methods is proved. However, there is a subtle point: it is not clear a priori 
whether a formal connection between static and TAP approaches is valid only for the correct full RSB solution of 
the SK model, or if such a connection is preserved at each finite (though approximated) level of replica symmetry 
breaking. Our result shows that this second scenario is the correct one. 

A further, highly nontrivial consistency check of the static and TAP approaches, comes from a comparison of 
the TAP complexity with the complexity computed by means of constrained thermodynamics, which results in the 
Legendre transform construction of ]2C| ]. Within this method the complexity is given by the Legendre transform of 
the thermodynamic free energy of the system, thus establishing a connection between TAP and static methods at a 
generic value of /. We will prove that, as in the p-spin spherical model, in the SK model the TAP complexity, once 
the BRST relations are considered, coincides with the complexity of pof . 

In Section II we introduce the BRST supersymmetry in the context of the TAP approach, and derive the BRST 
relations that will be used in the rest of the paper. In Section II we show how the BRST relations enforce the DDY 
proof of (||), without the need of any extra hypothesis. In Section IV we present the BM calculation of the SK 
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complexity, exploiting the BRST supersymmetry, and in Section V we show the connection between our results and 
the Legendre transform method for the computation of the complexity. Finally, we draw our conclusion in Section 
VI. 

II. THE BRST SUPERSYMMETRY IN THE TAP CONTEXT 

In this Section we show why the BRST supersymmetry is helpful in the context of the TAP approach. This is just 
a specific example of the application of the supersymmetric formalism to statistical mechanics, and in particular to 
the field of disordered systems |24|, [27|, Q . As we have seen, the TAP free energy Ftap is a function of the 

local magnetizations rrii, and the minima of this function are identified with the metastable states of the system (TAP 
states). In the rhs of equation (|J) we have a typical example of a sum over different TAP states, as often needed in 
this approach. In addition to the partition function Ztap, we may also compute the density of TAP states ([!]), 

M 

p(f) = Y / S[F T Ap(m a )~Nf), (5) 
which is the quantity one needs for the computation of the complexity £(/). In (Q), as in (§), TAP states arc labeled 



by a, and m a indicates the corresponding set of local magnetizations. More generally, in t 
have to deal with expressions of the form 



ic TAP context we always 



M 

R=Y,r[F T Ap(m a )} , (6) 

Q = l 

r being a generic function of Ftap- The quantity R can be written as, 

R=Y, [l[ dm * 5 ( m i - O r[F TAP {m)] = 5{diF TAP {m)) | det(d^F TAP (m))| r[F TAP {m)] , (7) 

a—l i i 

where the TAP states have been identified with solutions of the TAP equations diFxApim) = 0. In ([?]) the modulus 
of the determinant is quite hard to handle, and due to this it is disregarded in most supersymmetric calculations 
p9j . This approximation is a priori unjustified, since without the modulus we are weighting each TAP solution with 
the sign of the Hessian determinant, with the risk of uncontrolled cancellations. The situation may slightly improve 
if the function r(FxAp) is peaked only on low values of the free energy Ftap, as it happens in the calculation of 
Ztap at low temperatures, and of p(f) at low free energies. In these cases, we may hope that at low temperatures 
(free energies) only minima of Ftap will dominate the integral in (Q). Minima have positive-defined Hessian, and 
thus the modulus becomes redundant. Moreover, the very identification of TAP states with TAP solution is only 
sensible if these solutions are minima, rather than generic saddles. Yet, the validity of this scenario should be checked 
carefully, as it can be done for the p-spin spherical model, where the effect of disregarding the modulus is completely 
under control, due to the fact that the free energy distribution of minima and saddles is completely known |BQ]. In 
the SK model we do not have such a precise information about the TAP free energy saddle points, and therefore 
dropping the modulus is quite risky. However, our aim here is to show that all the classic TAP calculations (where 
the modulus was disregarded) are formally consistent with the standard static approach, thus justifying a posteriori 
this approximation. 

We can use an exponential representation for the delta function and the determinant, 



/ +ioa rlr- ( \ 

Ilf^ expnTxidiFTApim)] 



/+oo ^ 
: [ ilr, 'I i ', exp 
-oo 



2J ipii>jdidjF T Ap{m) 



(8) 



(9) 



where {ip,ip} are anti-commuting Grassmann variables. In this way we can write p6j, 

R= [ VmVxV^ViP eS(m,*,$,l>) ; ( 10 ) 
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where the action S is given by, 

S(m,x,j},ip) = y^^XidiF TA p(m) + ^ip. l 'ip : jdid j F TA p(m) + w[F TA p{m)} . (11) 

i ij 

In the case r is an ordinary function we have w = logr, while if r(FTAp) — S(Ftap) then w = uFtap, where u is an 
imaginary integration variable to implement the 5- function (see next Section) . The measures in (|Io| ) include the sum 
over the indices and the constant prefactors. A key property of action ( pd| ) is its invariance under a generalization 
of the Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry [[l5[ |l6| (see also f26||): if e is an infinitesimal Grassmann 
parameter, it is straightforward to verify that (O) is invariant under the following transformation, 



Srrii — eipi Sxi — ~ew' ipi S^i = —exi Sipi =0 => SS = . (12) 

This generalization of the standard BRST supersymmetry to the case where w ^ has been first introduced in |22| , in 
the context of a TAP calculation for the p-spin spherical model. The BRST invariance does not depend on the explicit 
form of the function FxApim), but simply on the formal structure of action ([TT]) . In particular, it is essential that w 
is a function of the magnetizations mt only through FTAp{m). The fact that SS = under the BRST supersymmetry 
implies that the average of any observable of the same variables performed with this action must be invariant too. 
For an observable 0, 

(0(A)) = J VK 0(A) e 5(A) , (13) 

with A = {m, x, if), Vli we have 

(0(A)) - (0(A - SA)) = (50(A)) = . (14) 

This property can be used to generate some useful Ward identities. A fruitful choice is = mfipi and = xfipi, 
whose variation gives, 

k(m k - 1 i) l ip t ) = -{m^Xi) (15) 

fcfcJ-V&Vi) = (x$Xi) . (16) 

In particular, the case k — 1 gives the BRST equations already used in |22j , 

(&Vi) = -(miXi) [BRST1] (17) 
(w'^i) = (x lXl ) [BRST2] . (18) 

The BRST relations are crucial. Thanks to them we will be able to prove in general that any TAP average must give 
the same result as the equivalent thermodynamic average, and more specifically to reduce the TAP calculation of the 
complexity to the standard replica calculation of the free energy 



III. THE DE DOMINICIS- YOUNG CALCULATION REVISITED 

In 1983, De Dominicis and Young |l2) (DDY) introduced the statistical weight (||) for the TAP states and explicitly 
proved that under a key hypothesis the average of the replicated TAP partition function is equal to the average of the 
replicated standard partition function, i.e. 

J VJP(J)Z n (J) = J VJ P(J)Z^ AP (J) , (19) 

which is the practical way to write equation ([}]). This result is very important, because it proves the identity of the 
TAP and static approaches before any Ansatz for the overlap matrix is done. The hypothesis that DDY invoked 
to justify their computation is that the complexity of the dominant TAP states at any temperature is zero. More 
precisely, we can write, 

Ztap = f df p(f) exp[-0Nf\ = I df exp[-/3N(f - TS(/))] - exp[-/?iV(/* - T£(/*))] , (20) 
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where the dominant free energy /* is solution of the equation, 
DDY assumed as a necessary condition that 

£(/*) = . (22) 



In the present Section we argue that, although relation (|22|) is in fact verified in the SK model, it is not a necessary 
condition for the self-consistency of the TAP approach and for its equivalence with the statics. 



A. The problem and its solution in the DDY formulation 

Let us briefly summarize the arguments of DDY leading to condition ( p2| ) (in order to keep the notation as clear as 
possible we treat m as a one-dimensional variable in this Section; things do not change with the full iV-dimensional 
representation). First, we introduce an auxiliary magnetic field h both in FxApijn) through the term —hm, and in 
the Hamiltonian H(o) with —ha. From the very definition of equilibrium magnetization and energy we trivially have, 

-R M = E CT e M-PH(a)} = (a)Z , (23) 



P dh 
_dZ_ 

8(3 



]T H(a) exp[-0H{a)] = (H) Z . (24) 



On the other hand, from (||) we get, 



-a g h = X ~ l_ exp[-/3F TAP (TO Q )] = X + (a) Z (25) 

- Y + y — exp[-l3F TA p{m a )\ = Y + (H) Z . (26) 



d(3 ^ 5/3 

a 

To obtain this result we used the following relation, 

dF T Ap(m a ) 



a eM-PH(a)} , (27) 



dh ^ 

where the sum over a is restricted to those configurations belonging to state a (an analogous relation is valid for 



dFxAp/dP). On the other hand, in equations (25)-(p6[) X and Y are the contributions coming, respectively, from the 
dependence on h and f3 of the total number of TAP states, JV(h, 0). DDY correctly noted that X — and Y — are 
necessary conditions for the consistency of the TAP approach. Moreover, by using ((2^), they noted that if £(/*) = 0, 
these extra contributions are indeed vanishing in the thermodynamic limit, and thus consistency is recovered. Thus, 
DDY assumed £(/*) = and accordingly imposed X = Y = in their computation. More precisely, after averaging 
over the disorder, the relations X = Y = become two equations for the auxiliary fields (eqs.(32) and (42) of fjjl), 
which DDY used to prove (||). Our point, however, is that the equations X = Y = are always satisfied, irrespective 
of the value of £(/*). 



B. Role of the supersymmetry 

We now show that the mathematical conditions on the auxiliary variables that DDY impose are nothing else than 
the BRST relations, which are valid in general, and contain no information on £(/*). Let us define, 

G ( m , h) = d (j3 f TAp) , r(m, h) = e -^AP(m,fc) . (28) 
dm 
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We consider relation (|23|) (a similar reasoning can be easily done also for relation Q24J)) and use the integral represen- 
tation of Section II for Ztap, to obtain, 

= [dmdxx e xG ^ |^ r(m, h) + f dm dx e xG ^ £ . (29) 
oh J dm J dm oh 

In this expression the first term is X , whereas the second term is just (a) Z. Integrating by parts we have, 

X = [ dm 6(G(m, h)) ^- ■ (30) 
J dm 

At this point note that for a generic function r(m, h) there is no need for X to be zero, in much the same way as 



action (11) is not BRST invariant if r is not a function of Ftap- If, however, r(/i, m) — r(FTAp(i r n, h)), as in the case 



under examination, we have, 

X = [dm 6(G{m, h)) G(m, h) ® r = . (31) 
J db T AP 

Thus, the extra term coming from the differentiation of Af(h, (3) with respect to h vanishes if the TAP states are 
weighted with a function of the TAP free energy. A similar conclusion can be easily drawn for differentiation with 
respect to j3, which leads to Y = 0. In particular, these relations hold in the case r(FTAp) — exp(— ^F^p), which is 
the one analyzed by DDY. 

The important fact is that from a mathematical point of view, X = , Y = are consequences of the particular 
form of action (fill ), and of its symmetry properties. Indeed, if we express ( p3| ) and ( p4[ ) in the supersymmetric 
representation (|13|), we get, 

X = (x) 

9G, 



Y 



dm 



— l -y^c k {m k x) + k{m k - 1 ^), (32) 
dp ^ 



with, 



_ 1 dG 

Ck ~ k\ dm k d(3 ■ [ ' 

A comparison between (|3^) and (|l5|) clearly shows that the relations X = and Y = are a direct consequence of the 
BRST symmetry. Thus, the relations imposed by DDY in the calculation of ZJ^ AP , and in their demonstration of the 
equivalence between TAP and static averages, are a consequence of the BRST supersymmetric form of the TAP action 
( pi] ) . These relations are therefore not the expression for the absence of an extensive complexity of the equilibrium 
states, i.e. they do not imply £(/*) = 0. For this reason, the result of DDY is more general than what originally 
thought. This is consistent with the fact that a similar formal connection between TAP and static approaches is valid 
in the p-spin spherical model, where, in a certain range of temperature, £(/*) > 0. 



C. Disregarding the modulus 

Before closing this Section, a word of caution on the modulus of the determinant |27|, |30|| . As already noted 
above, disregarding the modulus can be very risky, especially if the function t(Ftap) weighting the TAP states in (|^) 
is not peaked on low free energies, and thus minima. As an extreme illustration of this risk we consider the function 
t(Ftap) = 1. Clearly, from (||) we have, 

m . m * f = ^°- (34) 

However, if we use for R the integral representation (Q) and we disregard the modulus, we obtain, 

f = °- < 35 > 

since the BRST supersymmetry is trivially satisfied with t(Ftap) = 1- The problem here is that by disregarding the 
modulus, and with a flat weight r, we are summing over all stationary points of Ftap, each multiplied by the sign of 
the determinant. The Morse theorem states that this quantity must be a topological constant, only dependent on the 
manifold over which Ftap is define and on the boundary conditions on Ftap- Result ( p5| ) is therefore correct, but 
the quantity in this equation is not the same R(h) as in ([m|). This trivial example shows how important is to weight 
TAP solutions with a function peaked as much as possible on low free energies when the modulus is disregarded. 
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IV. THE BRAY-MOORE CALCULATION REVISITED 

In the present Section we will compute the annealed complexity of the TAP states for the SK model, following 
closely the classic calculation of Bray and Moore (BM) . Our new contribution will be to exploit the BRST relations 
( |l7| ) and ( |l8| ) in order to simplify the resulting saddle point equations. In this way we will prove that the complexity 
is intimately connected to the 1RSB static free energy. 



A. General definitions 



The TAP free energy for the SK model is given by 



F TAP (m) = ~ JijmiTnj + — ^2 <M<7, 



/??,■ 



(36) 



with, 



<f>a(q,m) = -(1 + m) log 



(1 + m) 



+ ^(1 - m ) l°g 



-(1 — to) 



1 f3 2 
= — log(l — to 2 ) + m tanh~ (to) — log 2 — (1 — q) 2 . 

The variables rrii are the local magnetizations, and q is the self-overlap of the TAP states, 



A-^'" : 



(37) 



(38) 



while the quenched couplings J are random variables with Gaussian distribution, 

P{J l3 ) = exp(-iV4/2) . 

The TAP equations and the Hessian of the free energy are respectively, 

(3diF TA p{m) = -j3 y j Jijirij + (j>\(q, mj) = , 

(3did 3 F T Ap{m) = -fMij + ^{q^rrii) % , 



with, 



4>i{q,m) = /3 2 (1 — q)m + tanh 1 (m) 

1 



1 — m 2 



0(1/N) 



(39) 



The term of order 1/N in 4>2(q> m) will be dropped in what follows. Following BM we perform an annealed calculation 
of the number of TAP states, i.e. we directly average p(f,j3\J) in (@) over the distribution of the quenched couplings 

Jij, 

Af(P,f)= [ VJP(J)p(f,(3\J) . 



(40) 



We use for p the integral representation of eqs. (tL0|) and (|l l|) , with 



S(F TAP -Nf) 



(hi 



2ni 



t exp[u(F TAP -7V/)] , 



and thus w(Ftap) = u (Ftap — Nf) and w' = u. Thus, the average number of TAP states becomes, 

Af(f3,f) = f VJP{J) VmVxV^V^du e ^im,x^,u) ^ (41) 
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where we have multiplied the action by f3 in order to keep our calculation as close as possible to BM. Before proceeding 
we mention an important point: in the calculation of BM the J-dependent part of FTApijn) m S is eliminated by 
using the equations diFTAp{m) = 0, which are enforced by the (5-function. More specifically, BM use the equation, 

ij i 

which is valid in the TAP states. This substitution simplifies considerably the calculation, but unfortunately the 
action obtained in this way is no longer BRST invariant. Thus, in the present calculation we must use the full form 
of FxAp( m ), equation (|3q), and for this reason some parts of the calculation differ from BM. 



B. The calculation 

The J-dependent part of the action in ( [fl| ) is given by, 

[3S.j = — /?y^ J l j{x i rrij + ^ipj + \-u mirrij) , (43) 



and after averaging J\f over the disorder we obtain the new effective action, 

(3S = y^x a 0i(q, TUj) + y^ j '4>iipi^2(q,m i ) + </> n (g,mj) - N/3uf 

i i i 

+ ^5>? + ^ £(w) 2 - t^Ew^) 2 + ^t uV + ^ £ miXi ■ (44) 

i i i i 

In order to linearize the quadratic terms, we introduce in ( p~T| ) the following 5-functions, 

6(qN -£m?) = / d\ e- XqN+x ^> m > (45) 



— too 
-\-ioo 



5(RN -^m^) = / dr e- rRN+r ^* m ' x ' (46) 
5(TN-J2^i) = / dt e- tTN+t ^-^ , (47) 



and we integrate over g, R and T. In this way the integrals in xi and (tpi, ?pi) become Gaussian and can be performed 
explicitly. The effective action becomes, 



/3S(Cl,m) =iVS (O) +££(fi,mi) , (48) 



where 51 = {q, A, r, i?, t, T, u}, and, 



So(fi) = ^R 2 -^T 2 -rR-tT-hog(2TTp 2 q)-uff3-\q+^u 2 q 2 + f3 2 uqR, 

C(n,m) = u<t>o{q,m) - — L- [^(^m) +rrn] 2 + log[0 2 (<j, m) + i] + Am 2 . (49) 
We can now write the average number of TAP states as, 

mpj) =J Vil e N *°W Y[ J dm t = j ' TXl e ^o(o )+ ivio g / dm (5Q) 

Thanks to the prefactor N in the exponential, the integral in T>fl can be performed with the steepest descent method. 
In this way we can write the complexity E(/3, /) as, 

E(A /) = ^ logATpJO = E (A) + log / ^ e c ^ , (51) 
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where is solution of the saddle point equations, 



on 



with, 



((0(m))) 



1 



dm 0{m) e c( - Q ' m ) 



J dm e £ ( n ' m ' 

We can reduce the number of variables by directly solving the saddle point equations for R and T, 



5So 
dR 

dT 



R = r/f3 2 -qu 



T 



-t/P 2 



In order to have expressions as similar as possible to the ones of BM, we define, 

B = (3 2 (l-q)+t 
A = -/? 2 (1 -q)-s. 

Using the explicit forms of 4>ii4>2 in (|39|), and relations (Q) and (|55|), we can rewrite ( ^9[ ) as, 

(B 2 - A 2 ) 



s (r>) 



-\q-(3uf-(B + A)(l-q) + 



2{3 2 



- 1 \og(2nf3 2 q) - ^u 2 q 2 - uqA - (3 2 uq(l - q) 



£(fi, m) = log 



1 



B — 



[tanh (m) — Am]' 



+ Am 2 + u(j) (q, to) , 



1 - to 2 / 2/3 2 g 
where </>o is given in (|37j). The remaining saddle point equations ( fi^ ) for the variable {A, u, B, A, 5} are, 

q=((m 2 )) 

P \q 2 -qA- (3 2 q(l - q) 



^=0 
dX 



du 







/?/ = «&,(<?, m))> 



as 

<9A 
<9g 



A 



l-/3 2 (( 



(1 - TO 2 ) 



2\2 



1 + 5(1- TO 2 







'Y^ 1 ~ ^ + ^(( mtanh_1 ( m ))) ~ Y uq 



A = fl + A- i|l- ^(([tanh- 1 (m) - Am] 2 )) } 

d<p [q,m 



-))- 



fl2 

^-M 2 g + uA + /3 2 u(l - 2g) 



(52) 



(53) 



(54) 



(55) 



(56) 
(57) 

(58) 
(59) 
(60) 
(61) 

(62) 



These equations, together with relations (q6|) and ( |57j ) can be compared with equations (15), (16) and (17) of BM 
. The differences are due to the different representation of Ftap (jn) we have taken, in order to preserve the BRST 
supersymmetry. However, the final result, i.e. the values of S and C in the solution of the saddle point equations, is 
exactly the same (in the last of equations (17) of BM, however, there is a term udf jdq missing). Of course, if we set 
u = 0, i.e. if we do not impose the constraint on the free energy, our expressions become formally identical to BM. 



C. Role of the supersymmetry and connection with the statics 

Solving the saddle point equations above, even numerically, is not a simple task. However, the problem becomes 
much easier if we make use of the two BRST relations jl7| ) and (|l8|). From (|47|), ( Eij ) and (pq), we have, 



(iPi^i) = T = ~t/(3 2 = -BI0 2 + (1 - q) 

{x l m i ) = R = r/0 2 - qu = -A//3 2 - (1 - q) - qu 
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and thus the first BRST equation becomes, 

B + A = ~(3 2 q u [BRST1] . (63) 

In order to write the second BRST equation we need {xiXi}, which can be computed from equation (pfj]), 

te> = ~^r q j 1 - ^((t tanh_1 ( m ) ~ Am l 2 » 

In this way the second BRST relation becomes (we recall that w' = u) , 

- I 1 " -4-(([ tanh_1 M- Am ] 2 ))\ =u[B- I3 2 {l-q)] [BRST2] . (64) 

q I p q J 

Equation ([30|) admits the solution B — 0, and following BM this is the solution we adopt. Setting B = into the first 
BRST relation we obtain, 

A = ~[3 2 qu , (65) 
while substituting equation ( |37| ) and the second BRST relation into (|62|), we find, 

A = i/?V<z . (66) 

We now use these two results, eqs. ( |65| ) and (|66|) , to rewrite the complexity E(/3, /) in (pT|). By making the change 
of variable m — ► ft, = tanh _1 (m) in the integral in (|5l|), and using (E^) and (B^), we obtain, 



E(/3, /) = tt |- log2 - i- log(2^/3 2 g ) + £-[(-« - l) g 2 + 2g- 1] + ~log^ dh e ^9>«) j - 



(67) 



with, 



ft 2 

T{h;q,u) — — — w log cosh ft,. (68) 

2p^q 

This form of the complexity can be fully appreciated by recalling the expression for the free energy of the SK model 
at the 1-step level of replica symmetry breaking (1RSB). The 1RSB free energy density is given by 



pF 1RSB (P;qi,x) = - log 2+ i-log(27r/3V) + ^[{x - l)ql + 2 9l - 1] - ~ log J dh e^^i,-*) , 



(69) 



where the self-overlap q\ and the breaking point x satisfy the saddle point equations dFi R $ B / dqi = 0, dF\ R s B /dx = 0, 
and where we have set to zero the mutual overlap qy. Clearly, we see that there is a striking formal correspondence 
between E and FmsB- Indeed we have, 



E(/3, /; q, u) = (3u [F 1RSB (ft q, -u) - /] , (70) 



and, 



n dYl a dF 1RSB (P;q,-u) 
= ^ = PU dq 

<9E dF 1RSB ((3;q 7 -u) 

~du = du \- P[FiRS B {p;q,-u) - /] . (72) 

Therefore the saddle point equations for E and Fi R $ B coincide, provided that / = Fi R s B (f3;qi,x). Moreover, for 
/ = F\ R s B we trivially have 

E[/3,F lfl ss(/3;?i^)] = 0, (73) 

and therefore /o = Fi R g B : as expected, and anticipated in the introduction, the lowest TAP states have the static 
free energy density. 
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The fact that Qi= q and x — —it, at fn can also be obtained by directly setting f = fo = F\rsb in the saddle 
point equations ([38||32|). From ( |59|) and (|6l]), after some algebra, we obtain, 

--/3 2 q 2 u + — log(2n(3 2 q) - - log / dh e ^ h ' q '^ - ((logcosh/j)) = , (74) 
4 2u u J 

while equation ([38|) becomes, 

l + ^qu[(u+l)q~l]--^((h 2 )) = , (75) 

where the averages ((•)) are now performed with the distribution exp[jF(ft,; q, u)] of (|68|). Once we set it = — x and q — q\ 
equations (|74|) and ( |75| ) coincide with the 1RSB saddle point equations, i.e. OFirsb /dx = and OFirsb / 'dqi = 
y. Summarizing, the complexity S vanishes at the 1RSB free energy density, and at this point the two calculations 
are formally equivalent, since the saddle point parameters coincide, and £ and Firsb are related by equation (|70|). 

Some considerations are in order. First, this identification of the saddle point parameters in the two calculations 
has a clear physical meaning, already discussed in p^ ]. The self-overlap qi of the 1RSB states is just the same as the 
self-overlap of the TAP states at / = Firsb- Less obvious is the relation u = —x, since the 1RSB breaking point x 
does not seem trivially related to the parameter u. However, we must note that, 

^ = -M/) , (76) 

and thus —Bu(Firsb) is the slope of the complexity at the lowest free energy, which is indeed equal to the static 
breaking point [3x. 

Secondly, we recall that we have set qo = in the 1RSB calculation. This is due to the fact that we performed 
an annealed computation of E, and thus we only had one value of the overlap. We believe that a quenched, but 
replica symmetric, calculation of T, will be equivalent to the 1RSB static calculation with q$ ^ 0. This brings us to 
the final point: the TAP method has one less step of replica symmetry breaking as compared to the standard static 
one. This is simply due to the fact that the elementary objects in the TAP approach are states, while in the static 
one are configurations. For this reason, a model as the p-spin spherical spin glass, which is exactly solved by a 1RSB 
Ansatz, has a complexity which is exactly replica symmetric. In the SK model, however, which needs a full RSB 
static solution, also the correct complexity will need to be computed at a full RSB level p2| . 



D. Numerical solution of the saddle point equations 

We now consider the saddle point equations for £ at a generic value of / > Firsb- We are left with two unknown 
variable, q and u, and with three unused equations, that is (p8[), ( ]59| ) and (|6l|). However, by using the BRST relations 
it is possible to show that actually equations (|5^) and (^TJ) coincide. To prove this fact one must use the formula, 

((h 2 )) = (3 2 q - (3 i q 2 u + (3 4 q 2 u(u + 1) ((tanh 2 h)) , (77) 

in (|o4|), and then substitute this second BRST relation into either (|58l), or (|6l"l). Thus, the two remaining saddle point 
equations for q and u, (|58|) and (|5S|), become, 

q = ((tanh 2 /i)) 

B 2 

8f + ((logcosh/i)) + ^-[ g 2 (l + 2u)-2g+l]+log2 = . (78) 

where the formula, 

~ log(l — m 2 ) = — log cosh h , 

was used. It is not difficult to solve numerically the two equations above, and this shows how drastic is the simplification 
of the calculation due to the use of the BRST supersymmetry. In Fig.l we plot £ as a function of / — Firsb{I3), 
at T — 0.7 T c (we recall that T c = 1). The system of equations ( |78| ) have solution only up to a threshold value 
/ = fth(P), where the complexity takes its maximum value T, t h- The threshold free energy turns out to be quite close 
to Fij?sb(/3). An expansion of the equations close to T c shows that / t /j = Fxrsb + 0(e 5 ), where e = T C — T, whereas 
q = 0(e) and u = 0(e). 
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FIG. 1: Complexity E as a function of the free energy density difference / — Fxrsb(P), for T = 0.7T C . 

An interesting comparison with the result of BM can be done by computing the total complexity of TAP solutions, 
Titot- BM obtain this quantity by removing the constraint given by the delta function on the free energy, that is by 
setting u = in the saddle point equations and in the complexity. The complexity S tot that BM find in this way is 
nontrivial (Fig.l of 0). If we set u = in our equation j67|), however, we obtain a trivial result, that is E = 0. This 
is consistent with our initial remarks on the modulus of the determinant: without the constraint on the free energy 
we are counting all the stationary points of the TAP free energy, each one weighted with the sign of the determinant. 
This quantity must be a topological invariant due to the Morse theorem, and thus the complexity must vanish. The 
fact that this happens in our calculation and not in the calculation of BM shows that the BRST symmetry selects 
the saddle point solution which preserves Morse's invariancc. 

In order to obtain £ tot in our case we simply note that at each temperature the dominant value of the complexity is 
given by its maximal value, and therefore E tot (T) = S t ^(T). Thus, the (annealed) average global number of solutions 
of the TAP equations is given by Aftot(T) ~ exp(iVE^(T)). The threshold values fth and T,th can be obtained by 
requiring that the system of equations ( |7q ) ceases to have solution, i.e. by imposing a marginality condition on 
their Hessian. The quantity T, t h(T) is plotted in Fig. 2: it turns out to be considerably smaller than the non-BRST 
symmetric value of . This might be one of the reasons why it is so difficult to find solutions of the TAP equations 
numerically. A numerical analysis of £ tot (T) shows that it goes to zero for e = T c — T — > as, 

E tot (T) = 0.01 e 6 + 0(e 7 ) , (79) 

We conclude this Section with an interesting remark on the comparison between the SK and the p-spin spherical 
model. In the p-spin model when we compute the complexity without the constraint on the free energy (i.e. setting 
u = 0) and without the modulus, we find two saddle points: a BRST symmetric one, which gives E = 0, as in the SK 
case, and a non-BRST symmetric saddle point, which gives a nontrivial complexity, E' 7^ 0. In the p-spin spherical 
model this second value E' of the complexity coincides with E^, that is the maximal value of E(/) obtained by 
keeping u 7^ and which is BRST symmetric. The fact that our result, Fig. 2, does not coincide with the result of 
BM (Fig.l of 0), shows that this identity does not hold in the SK model. Thus, the identity E' = £ t /j seems to be a 
peculiarity of the p-spin spherical model. 



E. The zero temperature limit of the total complexity 

The largest value of the total complexity £ tot (T) is achieved at T = 0. We find, 

lim E tot (T) = 0.0073 . (80) 

T-v0 

This value is remarkably smaller than the value E° ot = 0.199 obtained in 10, [5l] by directly solving the zero- 
temperature limit of the TAP equations. In this Section we want to discuss this difference and its implications. 



13 



0.008 



0.006 



W~ 0.004 



0.002 




FIG. 2: The threshold complexity T, t h as a function of the temperature. Due to the exponential form of the number of states, 
Etft is also the total TAP complexity of the system. 



The TAP equations at finite (3 are given by, 

mi = tanh{/3[/ii - (3(1 - q)rm]} 



where the local field hi is, 



JV 



For very large, but finite, (3, we can write, 



hi X! "' • 



= sign[/i, - (3(1 - q)rrii} , 



(81) 



(82) 



(83) 



where the magnetizations are now spin variable, m,; = ±1. In the limit (3 — ► oo we have q — » 1, and thus we must be 
careful with the expression in the square bracket. Let us write in general for (3 ^> 1, 



such that the TAP equations become, 



(1 - q) = h p~ 



rrii = sign 



a > , 



h mi 



(84) 



(85) 



According to the value of a we have different scenarios. For a < 1 there is no finite limit of the TAP equations 
This case is therefore uninteresting. For a = 1, on the other hand, we have a well defined zero temperature limit of 
the TAP equations, i.e 



= sign [hi - h rrii] 



(86) 



This is our case: in figure 3 we show that in the limit T — > the overlap q corresponding to the threshold (i.e. total) 
complexity behaves like, 



q = 1 - h T 



h = 0.49 , 



(87) 



so that a = 1. Thus, the zero-temperature limit of the BRST complexity, E" ot = 0.0073, is the complexity of the 
solutions of equations ( |S6| ) . Note that the T-dependence of the overlap may be different from ( |8~7j ) at values of the 
free energy different from the threshold. 
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In the case a > 1 the zero temperature limit of the TAP equations (|8^) yields, 

nii = sign [hi] . 

These are the equations considered in [l0| |TJ in order to compute the zero temperature complexity (in particular, 
equations (|88| ) were obtained by assuming that q = 1 — 0(T 2 ), i.e. a = 2). Therefore, the value E° ot = 0.0073 we 
obtain as zero-T limit of the BRST complexity and the value S° ot = 0.199 obtained in [HJ |3l]], refer to two different 
sets of equations, that is ( |S6|) and (^). The fact that solutions of ( |S6| ) have a much smaller complexity than those of 
(|S8|), is due to the fact that the former are actually a subset of the second. Indeed, from equations (Q) we have, 

\hi\ > ho , (89) 

for all i, which is a very restrictive condition not necessarily satisfied by the solutions of (|88|). 
We note that equations (|8^) imply, 

m ^ > ' ( 9 °) 

and thus, given a TAP solution, the change in energy AE(K) when we flip K spins, satisfies the relation, 

K K K ( h \ \ 

AE{K) = ^m i h i -^J ij m i m j ~Y^mihi-—== > K ( - —= J . (91) 

i ij i * \ ^ V / 

The condition of stability under K spin flips ]33| requires AE(K) > 0. A sufficient (but not necessary) condition for 
this to hold, is given by, 

k ° 1 >0. (92) 



In the case a — 1 this condition is trivially satisfied, while for a > 1 the order of the two limits (3 — > oo and N — * oo 
becomes important, and stability under K spin flips may be trickier to prove. Nevertheless, a word of caution on the 
concept of stability is in order at this point. At finite temperature a stable solution of the TAP equations is a local 
minimum of the TAP free energy, that is a solution with positive Hessian. Unstable saddle points of the TAP free 
energy will exist, typically at high free energies. Unfortunately, at T — this topological stability is impossible to 
assess, because of the discrete nature of the model. On the other hand, the -ftT-spin flips stability is easily defined at 
T = 0. However, it is not straightforward to understand what is the relation between these two definitions of stability. 
We simply note that we have computed the total complexity by finding the threshold states: for free energies above 
the threshold value the saddle points equations become unstable. Results in continuous mean-field spin-glasses (as the 
p-spin spherical model) suggest that the threshold may indeed be the border between topologically stable solutions, 
that is TAP minima, and unstable saddle points. It is therefore plausible that our total complexity refers only to 
TAP minima. The connections we have discovered between complexity and thermodynamics support this hypothesis. 
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V. CONNECTION WITH THE LEGENDRE TRANSFORM APPROACH 



As discussed in the previous Section there is a precise correspondence between the TAP computation at the lower 
band edge /q and the computation of the thermodynamic free energy of the system. However, the relation between the 
TAP approach and the usual static one (which involves Boltzmann averages) is still deeper: as it can be appreciated in 
© there is a formal connection between the TAP complexity at generic f and the free energy function Firsb ■ These 
equations also represent an important bridge between the TAP approach and some different methods to compute the 
complexity which do not rely on TAP equations, but rather on constrained thermodynamics p0| , pl[ . In pcfl it is 
argued that the complexity S(/3, /) of the ergodic components - i.e. states - present at low temperature is given by 
the Legendre Transform of the free energy F(f3, n) of n coupled real replicas of the original system. More precisely, 

E(/3, /) = max \fin f - /3F(j3, n)] , (93) 

71 

where n and / are Legendre conjugated variables, i.e. 

/ = ^M. ,94, 

If we assume, as it is generally accepted, that TAP minima correspond to low temperature states, then the complexity 
computed via eq.(|93]) and the complexity obtained from the TAP approach should be equal. In other words, we expect 
that the TAP complexity to satisfy eqs. (|93| ) and (04). 

To establish this point one has first to compute the free energy F(f3, n) of the coupled system with the replica 
method and choose an appropriate Ansatz. This problem is discussed for the SK model in [2q| , where different kinds 
of Ansatz are analyzed: within a generalized two-group Ansatz for F(f3, n), equation ( p3|) is formally satisfied by the 
BM solution, as already noted by these authors in j7|, g|. Moreover, a simpler sub-class of two-group saddle point 
solutions exists, which corresponds to the following procedure: F((3, n) at the fc-level of RSB is given by n times the 
free energy of one system evaluated with k + 1 steps of RSB, and breaking parameters nx\, nx2, ■ ■ ■ nxu, n. Thus, at 
the lowest possible level of RSB (A: = 0) one has, 

F((3,n) =nF 1RSB (P;q 1 ,q ,n) , (95) 

as also suggested in [ p0[ . With this expression of F(f3,n), eg. (|93] ) has been verified within the TAP approach for 
models which are exactly solved by 1RSB replica solution |2^, |3^. It has also been derived for the Bethe Lattice 
at zero temperature with the cavity method H ], and seems therefore quite robust. Surprisingly, for the SK model 
this point was still unclear, since the BM solution M , despite its formal correspondence to a more general two-group 
static solution, does not appear to satisfy relation (^3[ ) with the simpler Ansatz ([95]). The reason for this can be 
understood in the light of the previous Section: in tt BM do not use the BRST supersymmetry and consequently 
consider a larger set of solutions than those physically relevant, much as the generalized two-group Ansatz has a 
larger set of saddle points than the simpler Ansatz (^||). We therefore may expect that once the BRST relations are 
taken into account, as we have done in the previous Section, consistency must be recovered. Indeed, this is precisely 
what happens: as it can be easily verified, once (|9^) is implemented, eq.([7C|) is equal to ([93]), with qo = 0, q± = q and 
n = —u. The parameters q\ and n are variationally fixed: (|7^) is equal to (|94|), while variation with respect to q\ 
gives back equation (p7l|) . 

Thus, also for the SK model, we have demonstrated within the TAP approach the validity of relation (p3|). We note 
that a crucial element to establish this point is the exploitation of the BRST relations. This is not a surprise since, 
as discussed in Section III for equilibrium averages, the BRST supersymmetry is what mathematically guarantees 
relations of physical significance as (p3l) and (p3). 



VI. CONCLUSIONS 



In this paper we have discussed the equivalence between the TAP approach and the standard thermodynamic 
method for mean-field spin-glass systems, and we have shown that the use of the BRST symmetry of the TAP action 
is crucial to establish this equivalence. 

We demonstrated in Section III that the BRST relations guarantee the mathematical consistency of the TAP 
approach and ensure the validity of physically relevant relations. We have shown that the DDY computation, which 
establishes for the SK model a formal connection between the TAP and the thermodynamic partition functions, is 
generally valid if the BRST symmetry is taken into consideration, and does not rely on any hypothesis about the 
nature of the equilibrium states. Furthermore, we have revisited the BM calculation of the complexity for the SK 
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model, with the support of the BRST relations. Thanks to them, we have been able to solve the saddle point equations 
at given free energy, we have explicitly shown the equivalence of the lowest TAP solutions with the thermodynamic 
equilibrium states and we have verified the consistency of the TAP computation of the complexity with the Legendre 
Transform construction. 

It should be noted that the BRST relations have a very delicate role due to the presence of quenched disorder. 
Indeed, once the average over the disorder is taken and the auxiliary variables are introduced, the BRST symmetry of 
the effective action is much less transparent. In the thermodynamic limit, the auxiliary variables become variational 
parameters which must satisfy the saddle point equations. The solutions of these equations correspond to averaged 
quantities (for example, in the BM calculation R = (a^mi)) and should therefore satisfy the appropriate BRST 
relations. However, these relations are not automatically satisfied by all the saddle point solutions. In other words, 
the set of solutions of the saddle point equations in general includes also solutions which are not BRST invariant. 
Thus, the BRST relations must be explicitly imposed, providing some additional equations that must be satisfied by 
the variational parameters. Besides, as we have seen, this procedure provides a great simplification in the computation, 
since it drastically decreases the number of unknown variables. 

Our BRST calculation gives a total complexity at T — which is much smaller than the value previously obtained 
by other classic calculations. We have shown that technically this is due to the fact that the BRST saddle point 
provides a set of TAP equations at T = which are different from the ones previously considered. In particular, the 
zero-temperature BRST- TAP solutions have local fields larger than a finite value ho at all sites, and this condition is 
responsible for the drastic decrease in the complexity. We have outlined what are the possible physical explanations 
of this difference, but more study on the T — ► limit is necessary for a complete clarification of this point. 
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